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Abstract 

The distance from a given position toward one or more destina- 
tions, exits, and way points is a more or less important input variable 
in most models of pedestrian dynamics. Except for the special case 
when there are no obstacles in a concave scenario - i.e. each position 
is visible from any other - the calculation of these distances is a non- 
trivial task. This isn't that big a problem, as long as the model only 
demands the distances to be stored in a Static Floor Field also called 
Potential Field [TJ [2 El HI El El El [HJ, which never changes throughout 
the whole simulation. In this case a pre-calculation once before the 
simulation starts is sufficient. But if one wants to allow changes of 
the geometry during a simulation run - imagine doors or the blocking 
of a corridor due to some hazard - in the Distance Potential Field, 
calculation time matters strongly. This contribution gives an overview 
over existing and new exact and approximate methods to calculate a 
potential field, analytical investigations for their exactness, and tests 
of their computation speed. The advantages and drawbacks of the 
methods are discussed. 



1 Introduction 

1.1 General Introduction 

The will to move through space is the will to reach some kind of destination. 
On a smaller time scale this is the will to reduce the distance toward some 
kind of destination. Therefore it appears to be natural to use the gradient of 
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the distance toward a destination in some kind of measure as primary input 
and impetus for motion in the simulation of the movement of pedestrians. 

It is assumed that it is always sufficient to have a discrete distance poten- 
tial field, either because the model itself is formulated in a discrete manner, 
or because some finite - yet arbitrarily large - exactness of the distance po- 
tential field is sufficient. Although potential fields are by no means confined 
to rectangular grids [9], only rectangular grids are investigated. 



1.2 Robots and Pedestrians 

Note that the potential discussed in this contribution has the meaning and 
is used in the way as in the simulation of pedestrian dynamics [31 HI El 
EJ [10] and not robotics [TJ [2] . The difference is that in robotics it is usually 
assumed that an autonomous robot knows about his destination coordinate 
but has no knowledge of the position of obstacles except for those which it 
"sees". For pedestrians on the contrary it is typically assumed that they 
have at least some knowledge on the whole path, the positions of obstacles 
and the detours they have to walk compared to linear distance, even if they 
actually only see a small fraction of the whole path. This has consequences 
for the calculation and use of the potential. 

In robotics the artificial potential at position x from a destination at Xd 
was originally [TJ calculated as 

■tificial 

(x) = U 3d (x) + U obstades (x) (1) 
Ug d (x) = l -k d (x-x d ) 2 (2) 

leading to a potential as shown in figure [TJ In such a potential local minima 
can occur and a robot has to be equipped with the ability to realize it is in 
a minimum and how to get out of it. 



Figure 1: Example for Ug d {x) from equation ([2]). 

Much different from this is the basic assumption in many models of 
pedestrian motion that pedestrians have a good global knowledge of the 
exact Euclidean distances and know about the shortest path between their 
current position and their destination under consideration of all obstacles. 
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This assumption can be unrealistic for very complex geometries like huge 
mazes (compare figured), but for many situations it comes close to reality. 




Figure 2: For pedestrians the static floorfield or distance potential is cal- 
culated differently to account for the global knowledge pedestrians have in 
most situations. Mazes pose an exception inasmuch as the essence of their 
existence is that it is hard, if not impossible, to acquire perfect global knowl- 
edge. 

Note: the notion Distance Potential Field was chosen to clarify that it is 
nothing but a direct look-up table for the distance, whereas the Static Floor- 
field is inversely proportional to the distance and therefore - although edited 
in a trivial way - a quantity that's derived from distance to destination. 

2 Methods for the Calculation of a Distance Po- 
tential Field 

Short Mathematical Parenthesis: Vector Norms 

In two dimensions the so-called p-norms are defined as 

\\x\\ p := {\Ax\ p + \Ay\ p )p , with (3) 

of these p = 2 is the well known Euclidean Metric (Pythagorean Theorem) . 
So, for p = 2 the norm has the everyday life meaning of the word "dis- 
tance". At first two other well known norms, which are relevant here, will 
be discussed: the ones with p = 1 and the limit p — > oo. In principle the 
problem can be stated this way: only the metrics for p = 1 and p —* oo can 
be calculated by a flood fill, what one wants, however, is the metric with 
p = 2 for which no equally fast and simple calculation method exists. 

In this paper the following notation is used: the distance D from a certain 
point to the exit on the shortest path for a method X is called D x . It is 
composed of a sum of distances df of (straight) line segments of the visibility 
graph [TT|, i.e. they connect sequently the exit, corners of obstacles, and 
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the coordinate under consideration. Note that for the flood fill methods the 
visibility graph does not have to be known explicitly for the calculations. 

2.1 Flood Fill Methods 

In Flood Fill (sometimes also called Wavefront) methods the distance is 
calculated by sequently moving cell to closest neighbor cell and by that 
summing up the distances. 

2.1.1 Manhattan Metric 

For p = 1 one has the famous Manhattan Metric - also called Taxicab 
Metric or Manhattan Distance. It was introduced by Hermann Minkowski - 
but note that the name "Minkowski Metric" is reserved for the elementary 
metric of special relativity. Just as a taxi driver in Manhattan needs to 
sum up the number of Streets and the number of Avenues which he has to 
cross during the drive to get an estimation of the distance, the distance D M 
according to the Manhattan metric simply is the sum of the absolute values 
of the differences in x- and y-coordinate. Flood Fill therefore only acts 
within the von Neumann Neighborhood, i.e. grid cells need to be connected 
via a common edge (see figure |3(a) ) . 

Ax = \5xj\ and Ay = \ 6yj\ (4) 

i i 

D M = *£df = Ax + Ay (5) 

i 

2.1.2 Chessboard Metric 

For the limit p — > oo one arrives at the Maximum Norm - also called Cheby- 
shev Distance D c or Chessboard Metric. It was introduced by Pafnuty 
Chebyshev and got its alternative names from the way the king in chess 
is allowed to move, respectively the fact that the distance measured with 
this norm is the maximum of the differences in x- or y-coordinate. In other 
words: flood fill acts within the Moore Neighborhood and the added value 
is the same (typically 1) whether grid cells are connected via an edge or a 
corner (see figure [3(b) ) . 

df = max(\6xi\, |%|) (6) 
D C = £d? (7) 

i 

The advantage of Manhattan and Chessboard metric is the fact that the 
typically very fast flood fill procedure can be applied. For 2 < p < oo it is 
in general not possible to calculate the correct distance (i.e. the Euclidean 
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(a) Manhattan metric 



(b) Chessboard metric 



Figure 3: Manhattan and Chessboard metric. 



distance) between two arbitrary grid cells with a flood fill with local rules. 
The major drawback of course is that aside from Manhattan, Mannheim, or 
chess there are only few occasions where one of the two norms gives an exact 
result. However, there are possibilities to stick with the flood fill method 
and gain exactness. Some of these will be discussed in the following. 



2.1.3 Variant 1: Combination of Manhattan and Chessboard 

There is a method which uses flood fill and which is exact in the Euclidean 
sense for all positions which are visible from the destination. Since Manhat- 
tan metric gives df^ = |5xi=i| + |<%i=i| and Chessboard metric results in 
df =1 = max(|5xj=i|, |<%i=i|), one can simply calculate d 7 ^ =1 = dfL x — df =1 = 



mm(\8xi = i\,\5yi=i\). This "Minimum Norm" (see figure 4(a)) is of no use 

2 / \ ^ 2 2 

in itself, but (d^Li) + vH=\) m ust equal (5xi=i) + (<%i=i) , as one of the 
two 5s as trivially as necessarily needs to be the maximum and the other one 
the minimum. And (5xi = i) 2 + (5y i=1 ) 2 is the square of the exact Euclidean 
distance (p = 2) between the exit and a coordinate which is directly visible 
from the exit. If one generalizes this for line segments which are not directly 
visible from the exit (i > 1) one makes an error, since one can only calculate 
the square root of the sum of squared line segment (Euclidean) distances, 
while the exact result would be the sum of the square roots of squared line 
segment (Euclidean) distances: 



E( d f) 2 ^Ev( d f) 2 (8) 



Regardless of this problem, equations (|9J) (11) can be motivated by the 
initial observation that the calculation of distances is exact "to the point of 
the next corner" . This is achieved at the expense of making two flood fills 
- the final value of D Vl can not be calculated by a single flood fill - and 
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having to calculate a square root for each cell in the end. 

d? = df-df = min (\Sxi\, \S yi \) (9) 
D m = Y, d T = D M - D c (10) 

i 

(D v ^f = (D c ) 2 + (D m ) 2 (11) 



(a) "Minimum metric" (b) Variant 2 

Figure 4: "Minimum metric" and metric for variant 2 

2.1.4 Variant 2: V2 over Corners 

Another simple modification is a Chessboard flood fill where flooding across 
corners adds a \/2 instead of 1. If one does so, not only distances parallel to 
the discretization axis will be exact, but also distances deviating by exactly 
45° from that. A version of such a modification, which additionally includes 
a smoothing mechanism, is introduced in [12 . Since the number of diago- 
nal steps is min(5xi, 5yi) and the number of horizontal or vertical steps is 
max(5xi, 5yi) — min(5xi, dm) one can write 

df = V2dT + {df-dT) (12) 

D V2 = ^ d J 2 = (V2 - l) D m + D c (13) 

i 

2.1.5 Variant 3: Larger Neighborhoods 

A modification which gains computation speed on the cost of exactness is to 
increase the neighborhood further, so cells with a distance of two or three or 
even more are included. As one only has to continue the calculation with the 
border cells as new center cells, the recursion depth or stack size is reduced, 
but one runs the risk of overlooking small obstacles. This method is not 
investigated further here. 

2.2 Dijkstra's Algorithm on a Visibility Graph 

Another method is to try to find a subset of grid cells which form a Visibility 
Graph [TJ [HJ 03] (compare figure [5|. These grid cells are all grid cells 
which are necessary as navigation cell for at least one arbitrary grid cell if a 
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pedestrian wants to move from that arbitrary grid cell around the obstacles 
to the destination. Two nodes of the visibility graph are connected if and 
only if they are mutually visible. Once one has created such a visibility 
graph, one can calculate the distance toward the destination for all grid 
cells which are part of the visibility graph using Dijkstra's Algorithm |14j . 
After having done that one can calculate the distance from all other grid 
cells toward the destination by making use of the visibility graph and the 
distance information now contained within it. Note that strictly spoken the 
method used to measure the computation time is the one from [13] and not 
from [TJ 111 j . where the latter one is probably more efficient. 




Figure 5: Visibility graph of a simple geometry. 
2.3 Ray Casting 

Another very intuitive method is an iterated Ray Casting. In detail, the 
following algorithm is applied: 

1. Calculate the distance toward the destination for all grid cells which 
are visible from the destination. 

2. If possible, find that grid cell Xq, which is the closest to the destination 
of all those grid cells which have been assigned a distance, but which 
have at least one neighbor which is neither an obstacle, nor has been 
assigned a distance toward the destination, or has been assigned a 
distance, but a distance which is too large. 

3. Calculate the distance toward the destination of all cells which are 
visible from Xq. If there are cells which are visible to Xq and which 
have already been assigned a distance toward the destination, this 
distance is only overwritten, if the newly calculated distance is smaller. 

4. Repeat steps 2. and 3. until step 2. does not find any grid cell 
anymore. 

It is important for the calculation time how the check for visibility is 
done. Concerning the ray tracing part one can use the Bresenham Line 
Drawing Algorithm [15] . But it is important that this algorithm is not used 
to draw a line (cast a ray) each time one wants to calculate the visibility 
of a grid cell. It is better to draw a rectangle around the whole scenario 
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("border") and cast rays from the cell in focus to each of those border cells. 
All cells before the first obstacle are then marked as visible, all behind as 
"not visible" . Because of the discreteness, it might occur that the ray casting 
toward neighbored border cells gives different results for the visibility of some 
cell (when a cell is part of multiple casts). In that case, the cell needs to be 
marked as visible, because otherwise in scenarios with narrow spaces "blind 
spots" can appear, of which the algorithm would claim that they are not 
accessible at all. With this strategy the number of lines, one has to draw, 
only grows as the border size instead of as the area. 

2.4 Other Methods of Error Reduction 

One has to distinguish the distance error in the distance potential field from 
the speed error in a model. If the pedestrians move on a discrete lattice as 
well, the lattice structure can lead to errors in the speed, just as it leads 
to errors in the distance. The speed errors can be compensated for by 
making pedestrians suspend certain moves depending on the ratio of corner 
versus edge steps they did in recent moves. Such strategies are proposed for 
example in \12\ . 

3 Analytical Considerations 

3.1 Errors for Manhattan and Chessboard Metric 

The maximal errors of the two simple metrics are both: well known and 
trivial to calculate, nevertheless for the sake of completeness, they are given 
in the following. The absolute error compared to the Euclidean distance df 
of a single straight line element of some path using the Manhattan metric is 



with (fi being the angle between the connecting line of the two points and 
the x-axis. This simply sums up to a total error of 




(14) 
(15) 
(16) 



with df 



\J\5xi\ 2 + \5yi 




(17) 



which always lies between the boundaries 








(18) 
(19) 



with D E 




(exact total Euclidean distance) 



The maximal error arises from diagonal motion. 
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The corresponding values for the Chessboard metric are 

ef = df-df (20) 
= df (max (\ cos ifi\,\ sin (fi\) — 1) (21) 
E c = J^df (max(|cos<^|, |sin<^|) - 1) (22) 

i 

with the latter one always satisfying the relation 

(V(X5 - l)D E < E c < (23) 
Here as well the extreme value is reached for diagonal motion. 

3.2 Error for Variant 1 (Combination) 

To calculate the (extremal) errors, one has to deal with in variant 1, is a bit 
more complicated than it is for the two basic metrics above. The error is 

E Vi = D V! _ D E ( 24 ) 
^max(|(fai|, |%|)) + (^2mm(\5xi\,\5yi\j) (25) 

i i 

-Y^v\^\ 2 + \^\ 2 

i 

Remember that 

\8x{\ = df\ cos ((pi) | and \5yi\ = df\ sin ((fi) | (26) 

and use as abbreviation 

Mi = max (| cos (<pi) |, | sin (<pi) |) (27) 
mi = min(|cos(^j) |, |sin((^j) |) (28) 



in equation (25 ). 



EVl = jEE^f {M t M J + m l m J )-Y J df (29) 



« 3 



/E ( d f ) 2 + E E d f d f (M%Mj + m im] ) - df (30) 

i i jy^i i 

E d ? ) 2 + E E d f d f (MiMj + mi m 3 -1)-J2 dR31) 

i i i 



D 



1 . EiE 3 -#jrffrff (M^ +m^ -l) 

N (E^f) 2 



(32) 
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With the relations 



< mi < VOl < Mi < 1 
CL5 < minij + MiMj < 1 



one gets for the error 



D 1 



1 + V0.5 



(E^f) 2 



< £ yi < o 



(33) 
(34) 



(35) 



which is extremal for direction changes from horizontal/vertical to diagonal 
or vice versa. No error at all arises, when the new direction can be generated 
from the old by a reflection at one of the axis or diagonals. The left part 
will take the most extreme value, if all of the df are equal. With TV" as the 
number of single straight line elements one gets 



D 



1 + V0.5 - 1 



N-l 

N 



1 < E Vl < 



(36) 



This confirms the initial observation that there is no error for N = 1. How- 
ever, the error in this variant depends on the number of line elements and 
therefore the number of obstacles in - respectively the complexity of - the 
scenario. In the limit N — > oo the error can in the worst case (denoted by 
the hat) be 

E Vl = D E (0.5 - 25 - 1) » -0.159^ (37) 

3.3 Error for Variant 2 (a/2 over Corners) 

For variant 2 the error is 



E v i = ( V2D m + D c - D r 



a/2-1 J D rn + D 



c 



D 



D 



(38) 
(39) 



J2 ((v^- l) mmQSxillSyi]) + max(|(^|, \S yi \) - df) (40) 



with equations (27|) and (|28|) this is 

E v 2 



i 



v^-lU + Mi-l 



V2- 1] mi + a/1 



7717 



(41) 
(42) 



bearing in mind that < rm < \/0.5 one finds the maximum of each sum- 
mand at 



TTli 
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which corresponds to an angle of exactly <j> = 7r/8 = 22.5° and any corre- 
sponding angle in the other seven octants. The maximum error then is 



E V2 = - 2\/2- I J D E (44) 

« 0.082Z) B (45) 

which is slightly better than the maximum error E\ x « 0.099 for variant 
1 for N = 2 With this method there can never be a negative error. If one 
wants the error to vanish on the average of all angles, one can add the values 

Shv = ~ = —r^ r ~ 0.948 (46) 

a 8(V2-1) 

s d = v^- = V2—-H « 1.341 (47) 

a 8(^2-1) 

when flooding horizontally or vertically respectively diagonally. I.e. there 
is a global factor of ~ 0.948 multiplied to any distance, as Sd/shv = In 
this case the directions tpx z (of the first octant) with exact measurements 
are at 



— ±J—r 

vr A/ 4-2x72 vr 2 



sin ip 1>2 = ± x I ^= 2 ( 48 ) 



which is approximately (p\ ~ 9.55° and <p 2 ~ 35.45°. In the range between 
these two angles distances are measured too large, outside of this range too 
small. 



4 Computation Times 

4.1 Geometries 

The test geometries were - as shown in figure [6]- a "typical" room, a maze 
with 200 x 200 grid cells, a circle shaped room with a diameter of 996 grid 
cells , a square shaped room with a side length of 3998 grid cells, a room 
with a large column in the middle, and a ring, the latter two with the same 
size as the circle. 

4.2 Results 

The following tables give the results for computation time (standard PC) and 
the largest distance from the exit that was found by the method. Among 
the flood fill based methods variant 2 is different from all other flood fill 
methods regarding the fact that taking a square root for each grid cell is a 
necessary part of the calculation. For square roots and other elements of 
the calculation the calculation time not only depends on the algorithm but 
also on the details of the implementation. 
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Figure 6: Test geometries. Walkable areas are colored black, walls white. 
Note that the scenarios were largely different in size and just scaled differ- 
ently to fit the page (compare description in the text). 

In the visibility graph method, the visibilities were calculated via the 
area method. This led to the comparatively large calculation times. Though 
most of the time they were smaller than those of the ray casting method 
with visibility calculation by the area method. The only exception was the 
circle scenario. The reason for this is that the node generating method 
presented in [13] generates many unnecessary nodes, as the bending of the 
circle's border was so small that many local neighborhoods might also have 
been part of a convex corner, although the inner border of a circle is entirely 
concave. If the circle had been given as geometric object this could have 
been avoided. 

The ray casting method which uses the border method for visibility cal- 
culations was pleasingly fast, although never as fast as the Manhattan flood 
fill, and probably not yet fast enough to be applied each time step for many 
destinations. 

5 Conclusions 

In this contribution various methods for the calculation of distances in an 
obstacle filled space were investigated for their deviation from the true Eu- 
clidean distance and the time consumption for their calculation. Starting 
with the two well-known metrics - Manhattan and Chessboard - variants 
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Geometry: Typical 


Calc. 


Maximal 


Maze 


Calc. 


Maximal 


Method 


Time [si 

L J 


Distance 




Time [si 

L J 


Distance 


Manhattan 


0.04 


1246.00 




0.00 


7507.00 


Chessboard 


0.06 


960.00 




0.01 


7198.00 


Variant 1 (Comb.) 


0.12 


1001.70 




0.02 


7204.63 


Variant 2 [y/2) 


0.08 


1091.35 




0.01 


7325.97 


Ray Casting (edge) 


1.18 


1052.85 




72.72 


7361.20 


Visibility Graph 


11.00 


1052.16 




386.2 


7325.97 


Ray Casting (area) 


51.00 


1053.72 




1156 


7507.00 


Geometry: Circle 






Square 






Manhattan 


0.06 


1202.00 




2.44 


7994.00 


Chessboard 


0.10 


996.00 




2.69 


3997.00 


Variant 1 (Comb.) 


0.18 


996.24 




5.52 


5652.61 


Variant 2 (y/2) 


0.10 


1037.52 




2.68 


5652.55 


Ray Casting (edge) 


0.34 


996.24 




4.29 


5652.61 


Visibility Graph 


100.80 


996.24 




628.00 


5652.61 


Ray Casting (area) 


19.00 


996.24 




1606.00 


5652.61 


Geometry: Column 






Ring 






Manhattan 


0.02 


1993.00 




0.01 


1971.00 


Chessboard 


0.03 


1408.00 




0.01 


1386.00 


Variant 1 (Comb.) 


0.06 


1524.69 




0.02 


1504.40 


Variant 2 (y/Z) 


0.03 


1650.30 




0.01 


1628.30 


Ray Casting (edge) 


8.89 


1565.08 




7.36 


1542.24 


Visibility Graph 


47.80 


1564.94 




7.00 


1542.09 


Ray Casting (area) 


163.00 


1565.26 




22.00 


1542.39 



Table 1: Computation time and maximal distance for the geometries of 
figure [6j 
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of such flood fill methods were reviewed respectively introduced. Following 
that a visibility graph and a ray casting method were discussed. It was found 
that compared to the standard metrics it is possible to significantly reduce 
the error while sticking to the flood fill method by making subtle changes 
to or using combinations of the standard metrics. The errors that remain 
balance with the fact that their calculation is significantly quicker than the 
calculation of the two other methods which are - in principle - error-free. 
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